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Abstract. It is shown that the collective nonlinear interactions between intense 
neutrino or anti-neutrino fluxes and a dense neutrino plasma are governed by a multi- 
dimensional coupled cubic Schrodinger equation in which the interaction potential is 
positive or negative depending on the neutrino type. The cubic Schrodinger equation 
describes the splitting and focusing of intense neutrino beams due to the nonlinear 
excitations associated with the modifications of the individual neutrino energies in a 
dense neutrino background. 
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Neutrinos and anti-neutrinos are produced in the Sun and in supernovae during violent 
processes. The propagation of intense neutrinos in ionised media (dense plasmas) is a 
topic of current interest in view of the fact that neutrinos may be massive £Q . Because of 
their weak interaction with charged and neutral particles, neutrinos (or anti-neutrinos) 
can travel great distances without being affected appreciably by material obstacles. 
Massive neutrinos can account for neutrino oscillations and associated flavour change in 
solar plasmas |2HSj, while in supernovae and cosmology they may play a decisive role in 
reviving shocks for supernovae explosions, as well as in accounting for the missing dark 
matter in our Universe jl]. 

In the present paper, we employ a semi-classical approach to consider the nonlinear 
interaction between neutrinos and anti-neutrinos in neutrino plasmas. Using an effective 
field theory approach, a system of coupled nonlinear Schrodinger equations is derived, 
and it is shown that these equations admit splitting and focusing of intense neutrino (or 
anti-neutrino) beams due to the nonlinearities associated with the weak nuclear forces 
of the neutrinos/anti-neutrinos that are interacting with the background neutrinos. 
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Suppose a single neutrino (or anti-neutrino) moves in a neutrino-anti-neutrino 
admixture. The energy E of the neutrino (anti-neutrino) is then given by (see, e.g. 

mm) 



E = ^p 2 c 2 + m 2 c 4 + V±(r, t), (1) 

where p is the neutrino (anti-neutrino) momentum, c is the speed of light in vacuum, and 
m is the neutrino mass. The effective potential for a neutrino moving on a background 
of it's own flavor is 01101 ( see also 0HU) V±(r,t) = ±2^/2G F (n — n), while the 
potential for a neutrino moving on a background of a different flavour takes the form 
V±(r,t) = ±V2G F (n - n). Here G F /(hcf « 1.2 x KT 5 GeV~ 2 , where G F is the Fermi 
constant, n (n) is the density of the background neutrinos (anti-neutrinos), and the 
+ (— ) is for the propagating neutrino (anti-neutrino). From Eq. (JTJ we can, using 
standard techniques [12], obtain a Schrodinger equation for a neutrino (anti-neutrino) 
wave packet \P(r,t) JS] (see also Ref. ^1] for a similar treatment of neutrino-electron 
interactions) 

f + ^-^ = 0, ( 2 , 

where m e fr = rn/y, 7 = (1 — v 2 /c 2 )~ 1 ^ 2 , and f is the magnitude of the group velocity of 
the neutrino wave packet. 

In the case of self-interacting neutrinos and anti-neutrinos, moving with the same 
group velocity, we have n = (l^+l 2 ) and n = (l^-l 2 ), where and \P- are the neutrino 
and anti-neutrino wave functions, respectively, and the angular bracket denotes the 
ensemble average. Thus, in this case Eq. (J2J) turns into the coupled system 

*lt + f VV± T m * +{2) " (l ^ |2)) ^ = °' (3) 
where a = h/m e g and (3 = 2\/2G F /h for neutrinos of the same flavor. 

Suppose now that the neutrino pulses are coherent and uni-directional. If one 
species, the pump species (denoted by the index p), dominates the number density as 
compared to the signal species (denoted by the index s), i.e. n p 3> n s , we obtain the 
equation 

Ms ad 2 V s 
l lti + 2lh 2 

for the signal wave function where the driving nonlinear term is determined from 



M n ad 2 V„ 



It is well known that Eq. (j4^jl have modulationally stable dark soliton solutions ^H] , of 
which the fundamental soliton is 

% = !^ tanh (^J(3/a%z^ exp(-i(3^t). (5) 
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Figure 1. Contour plot of ln[n s (t, z)/n s o] (b = 1) for the splitting case. 
2, 




Figure 2. Contour plot of ln[n s (t, z)/n s o] (b — 1) for the focusing case. 



Following Ref. [E], we perform a perturbation analysis of Eq. (|4ajl for a given 
solution \P p (t,z) to Eq. (J4pj) . Separating \P s (t,z) into its real amplitude y/n s (t, x) and 
phase 9 s (t, z), we have, to lowest order in t, 

9 a (t,z)Ktp[ \V p (t',z)\ 2 dt'np\V p (0,z)\% (6) 
Jo 

where the last equality follows if the interaction takes place over a time-scale short 
compared to the time over which the pump amplitude changes. The signal number 
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density n s (t, z) is then 

r a3t 2 /d 2 nJ0, z) dnJO.z) d , , ,\i 
n,{t, z) « n s (0, z) exp [ - ( £ 2 ' J + ^ lnn s (0, z)J J , (7) 

where n p (t,z) = \\P p (t,z)\ 2 is the pump number density. Using the fundamental soliton 
(jSJ as the pump pulse, with \Pq = n p (0, 0) = n p0 , and a Gaussian distribution for the 
initial signal pulse, 

n s (0,z) = n s0 exp (-z 2 /a 2 ) , (8) 

where a is the width and n s o = n s (0, 0), we find that 

n s (t,z) = n s0 exp | - b 2 ( 2 - (3 sech 2 C - 2) sech 2 C - 2b 2 tanh C sech 2 C r 2 1 , (9) 

where £ = z/z , r = at/ z^, b = z /a, and z = (/3n p0 /a) 1//2 . We note that b, t and ( are 
dimensionless. The solution Q is plotted in Fig. ^ where the splitting of the collective 
neutrino state can clearly be seen. 

In the case of neutrino-neutrino or anti-neutrino-anti-neutrino interaction, Eq. (J3J) 
will be modified according to 

+ ~ mn2) + ^ 2 ^ = °' ( 10 ) 

where i = 1,2. For a dark soliton pump beam (j3J) and a Gaussian signal beam (JHJ), 
the approximate dynamics, analogous to the neutrino-anti-neutrino interaction, can be 
obtained from Q by letting r 2 -> -r 2 . The result is seen in Fig. El where the focusing 
is depicted. 

The above discussion is based on the approximate solution Q . However, numerical 
studies of the system (|4"qj) and (j4j|) show the same behaviour, when (jSJ) and (JH)) are used 
as initial conditions. Thus, a more complete solution supports the notion that the 
approximate analytical expression Q gives an accurate description of the dynamics of 
interacting neutrinos and anti-neutrinos. 

In conclusion, we have considered the nonlinear interactions between neutrinos and 
anti-neutrinos, mediated by a dense neutrino background of the same flavour. It is 
shown that the self-interaction between neutrinos and anti-neutrinos comes about due 
to their coupling via the weak nuclear force. Such interactions are governed by coupled 
cubic Schrodinger equations for the slowly varying envelopes of the neutrino and anti- 
neutrino wave functions. Analytical and numerical analysis of the multi-dimensional 
cubic Schrodinger equations reveal that a collective neutrino state can be split or focused 
due to the presence of a dense soliton background. The present results suggest that 
nonlinear excitations in neutrino plasmas can significantly affect the propagation of 
intense neutrino and anti-neutrino beams. The neutrino beam splitting can account for 
the energy loss of neutrinos, while the self-focusing renders neutrino burst propagation 
over a long distance. 
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